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Abstract 

We study the phantom-like behavior of a DGP-inspired braneworld scenario where 
curvature correction on the brane is taken into account. We include a possible mod- 
ification of the induced gravity on the brane by incorporating higher order curvature 
terms of Gauss-Bonnet type. We investigate the cosmological implications of the model 
and we show that the normal branch of the scenario self-accelerates in this modified 
scenario without introducing any dark energy component. Also, a phantom-like be- 
havior can be realized in this model without introducing any phantom field that suffers 
from serious difficulties such as violation of the null energy condition. 
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1 Introduction 



Observational data indicates that our universe is currently under accelerating expansion [1]. 
It seems that some unknown energy components ( dark energy) with negative pressure are 
responsible for this late-time acceleration [2] . However, understanding the nature of dark en- 
ergy is one of the fundamental problems of modern theoretical cosmology [3]. An alternative 
approach to accommodate dark energy is modifying the general theory of relativity on large 
scales. Among these theories, scalar-tensor theories [4], f{R) gravity [5], DGP braneworld 
gravity [6] and string-inspired theories [7] are studied extensively. One interesting model 
of modified gravity is the Dvah-Gabadadze-Porrati (DGP) braneworld model, in which our 
universe is a (3-|- l)-dimensional brane embedded in an infinite Minkowski bulk. This model 
has a very appealing phenomenology because it predicts that 4D Newtonian gravity on the 
brane is recovered at distances shorter than a given crossover distance Vc, whereas at scale 
higher than Vc the 5D effects become important due to leakage of gravity into the bulk. 
On the other hand, this model can explain the late-time acceleration of the universe in its 
self-accelerating branch without need to introducing a dark energy term [8] . Unfortunately, 
this branch suffers from some instabilities such as appearance of ghost degrees of freedom 
[9]. 

An alternative class of modified gravity models is the family of the string-inspired gravities 
by considering additional curvature invariant terms such as the Gauss-Bonnet (GB) term 
[10]. GB term arises naturally as the leading order of the a' expansion of the heterotic string 
theory where, a' is the inverse string tension [11]. Nojiri et al [12] showed that a particular 
dark energy solution can be obtained from scalar-Gauss-Bonnct cosmology. There is another 
version of the Gauss-Bonnet gravity namely the modified GB or F(G) theory [13] that can 
also play the role of gravitational dark energy. 

One of the essential tools for understanding the nature of dark energy is to detect the 
evolution of its equation of state parameter u = P/p, where P and p are pressure and 
energy density of dark energy component respectively. Recent analysis on the data from 
the Supernova, cosmic microwave background (CMB) and large scale structure (LSS) show 
that the cosmological constant fits well to the data [14], however current data also mildly 
favor an evolving dark energy with an equation of state u larger than —1 in the past and 
less than —1 today [15] with evolution across uo = —1 line in the intermediate rcdshift. If 
such a result holds on with a large number of the observational data, this would be a great 
challenge to the current models of dark energy. Firstly, the cosmological constant [16] as a 
candidate for dark energy will be excluded and dark energy must be dynamical ( note also 
that the cosmological constant needs a huge amount of fine-tuning) . Secondly, the case with 
u less than —1 which is often dubbed as phantom dark energy [17] introduces new theoretical 
facilities and challenges in this field. Phantom fields are a sort of scalar fields with negative 
sign for the kinetic energy term. In fact, phantom fields suffer from instabilities due to 
violation of the null energy condition, and a phantom universe eventually ends up with a 
Big Rip singularity [18]. Thus it seems natural to seek alternative approaches to realize a 



2 



phantom-like behavior without introducing any phantom field in the model. Phantom-like 
behavior is the growth of the effective dark energy density with cosmic time and in the same 
time, the effective equation of state parameter should stay always less than — f. In this 
regard, it has been shown that the normal, non-self-accelerating branch of the DGP scenario 
has the potential to explain phantom-like behavior without introducing any phantom fields 
on the brane ( for a number of attempts in this direction see [19]). The phantom-hke behavior 
of 4-dimensional f{R) gravity is studied in [20]. 

With these preliminaries, the purpose of this paper is to construct a class of DGP-inspired 
braneworld scenarios that curvature corrections are taken into account via incorporation of 
a Gauss-Bonnet type invariant term in the brane part of the action. Actually, the Gauss- 
Bonnet invariant should be considered in the bulk action. However, our goal here is to 
consider possible modification of the induced gravity on the brane by incorporation of higher 
order curvature terms such as the Gauss-Bonnet type terms. In fact, since Gauss-Bonnet 
term is an invariant in 4-dimensions, it takes contribution in the field equations if there is 
a coupling between scalar degrees of freedom on the brane and this invariant term. This is 
the main reason for incorporation of the field in /(G, 0). In the previous works done by 
other authors in this subject, generally an effective phantom phase is realized by adopting a 
phantom ansatz for scale factor. Here as we will show, the phantom phase is realized on the 
normal branch of the DGP-inspired model with curvature correction on the brane. We note 
also that generally one can show, by re-construction method, that f{R) gravity is equivalent 
to general relativity plus the scalar field [20]. We investigate the effect of higher order 
induced curvature correction in two alternative approaches: Scalar- Gaiiss-Bonnet gravity 
with F{G, 0) term in the brane action and modified Gauss-Bonnet gravity with f{G) term 
that G is the Gauss-Bonnet invariant. 

The paper is organized as follows: In section 2 we start from the action of the DGP- 
inspired F{G, (f)) scenario and derive the corresponding field equations. Then, we study 
the cosmological implications of this model. We show that in the presence of the curvature 
correction due to a Scalar-Gauss-Bonnet term in the action, the normal branch of the scenario 
has a self-accelerating behavior. The expansion history of this model will be studied and 
a comparison between luminosity distances in this model, the ADGP and ACDM will be 
performed. On the other hand, in this scenario we obtain a phantom-like behavior without 
introducing a phantom field neither in the bulk nor on the brane. In section 3 we consider 
a modified Gauss-Bonnet gravity in the DGP braneworld and we obtain the cosmological 
equations of the scenario. We show that this scenario can also realize a phantom-like behavior 
by adopting viable ansatz for the scale factor in the appropriate subspaces of the model 
parameter space. Finally, our summery and conclusions are presented in section 4. 



3 



2 A DGP-Inspired F{G, cp) Scenario 



We start with the action of a DGP-inspired F{G, 4>) scenario as follows 



5= / d'X^f^g^n+ [ d'x^g 

J bulk ^ J brane 



m 



^R + mlK + F{G,(j))+C, 



(1) 



Here with ^ = 0, 1, 2, 3, 5 are coordinates in the bulk, while with = 0,1,2,3 are 
induced coordinates on the brane. m| is the 5-dimensional Planck mass and TZ is the 5- 
dimensional Ricci scalar. Also, ml is the 4-dimensional Planck mass and K is the trace of 
the mean extrinsic curvature on the brane in the higher dimensional bulk, corresponding to 
the York-Gibbons-Hawking boundary term [21] defined as 



K 



1X1/ 



- lim 

2 e^+O 



+ 



y=-e 



K, 



y=+e 



(2) 



y is the coordinate of the extra dimension and the brane is located at |/ = 0. is the induced 
Ricci scalar on the brane, and the Scalar- Gauss- Bonnet term is defined as 



F{G, 4>) = --d^m - Vi<P) + f{<P)G{R). 
By definition, the Gauss-Bonnet invariant G{R) is given by 

G{R) — B? — ARfj^i^R'^" -\- R^i^aisR'^'^"'^ , 



(3) 



where all terms on the right hand side are functions of the induced curvature on the brane. 
In the action (1), Cm is the Lagrangian of the other matters localized on the brane. It is 
important to note that in the matter action, the matter is minimally coupled to the metric 
and not to the scalar field, making the Gauss-Bonnet gravity a metric theory which leads to 
conservation of matter on the brane. Note also that we have incorporated possible modifica- 
tion of the induced gravity on the brane by inclusion of the higher order induced curvature 
terms via the Gauss-Bonnet invariant. The corresponding energy-momentum tensor T^^^ is 
defined as 

2 5{^CJ 



rp{m) _ 



(4) 



Variation of the action (1) with respect to the scalar field, gives the equation of motion for 
the scalar field on the brane. 



VV-V^'(0) + /'(0)G(i?) = O, 



(5) 



where the prime denotes a derivative with respect to the scalar field 0. On the other hand, 
the bulk-brane Einstein's equations calculated from the action (1) are given by 



T^AB - \9ABn ) + miSA^SB'' 



R 



1 



1 



,u - -9,uR -g'"'f{cl>)G{R) + 2f{ct>)RR^^ 
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pa 



where 



S{y) = SA^SB^'r^Jiy), (6) 

(7) 



The energy-momentum tensor corresponding to the scalar field and curvature are defined as 
[12] 



and 



vp 



respectively. Equation (6) can be rewritten as follows 



where T^i, is the total energy-momentum on the brane defined as 



(10) 



T,,^--g^''f{(P)G{R)+2f{(P)RR^-'+4f{(P)R%R''f'+2f{(P)R^'^--R^^ 

(11) 
(12) 

(13) 



From equation (10) we find 



and 



Gab — T^AB — ^SabT^ — 0, 



Gfxi/ — Rpi/ pi' R 



for bulk and brane respectively. The corresponding junction conditions relating the extrinsic 
curvature to the energy-momentum tensor of the brane, have the following form 



lim 





y=+e 1 











Rpu — -^Qpug'^^Rap 



y=0 



(14) 



Using the 5D Codacci equation, the field equation of the bulk and the junction condition at 
the brane lead to the conservation of the total energy momentum tensor of the brane so that 



v^r^. = 0. 



(15) 



Note that, the matter sector of the energy momentum tensor T^^) can be described by a 
perfect fluid with energy density p^"^^ and pressure p("*) and satisfies the continuity equation 
by virtue of the Bianchi identity 



p("^) + 3i/(l+u;("^))p("^) =0. 



(16) 
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2.1 Cosmological Implications of the model 



To study the cosmology of a homogenous and isotropic brane in this setup, we consider the 
following line element 

ds'^ = g^ydx^'dx" + h'^{y, t)dy'^ = -n'^{y, t)dt^ + a^{y, t)'^ijdx'dx^ + h^{y, t)dy'^. (17) 

In this relation, is a maximally symmetric 3-dimensional metric defined as = 5ij + 



k 



where k = —1,0,1 corresponding to possible spatial geometries of the brane and 
= XiX'''. Now, choosing a Gaussian normal coordinate system so that b'^{y,t) = 1, leads 
to the following form of the junction condition (14) on the brane with the non- vanishing 
components of the Einstein's tensor in the bulk 



lim 

and 



dyU 



y=+e 
y=-e 



it) 



2nml 
ml 



na 



k 1 



lim 

e-»+0 



dytt 



2n^a^ ri^a 2a? 

ki 



+ 



n 



y=0' 



it) 



ml 
ml 



n^a 



+ 



2ml 



y=o 



(19) 



With these equations, the generalized Priedmann equation for the cosmological dynamics on 
the brane can be written as follows 



ml 



a^ 



{GB) 



(20) 



where the last term in the right hand side of this equation represents the dark radiation term 
and C is an integration constant. It is important to note that this term can be neglected 
at present time due to its fast decaying behavior. The energy density corresponding to the 
Gauss- Bonnet term is defined as 



(21) 



and the corresponding pressure is defined as 

p{GB) 



1^2 _ + 8^(i/V) + 16i/^0/'(0), 



(22) 



where a dot refers to derivative with respect to the cosmic time of the brane. The modified 
Friedmann equation (20) can be rewritten as follows 



k + p^^^^ 1 



3ml 



1 ^.pM+p(GB) 



3ml 



)■ 



(23) 



where £ = ±1 is associated to two possible branches of the scenario corresponding to two 



different embeddings of the brane in the bulk and Tc 



2m; 



^ is the crossover scale. In the 
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distance scale lower than this scale, gravity behaves as usual 4D scenario but in the distance 
scales higher than the crossover scale, gravity leaks to the extra dimension and this leakage 
leads to weakness of gravity in the large scales, so the universe accelerates. In other words, 
weakness of gravity at large scales plays the role of a negative pressure fluids that accelerates 
the universe expansion. Equation (23) as the modified Priedmann equation in this DGP- 
inspired F{G, 0) scenario is similar to the standard induced gravity result, but we note that 
here there is an important difference which is hidden in the definition of the total energy 
density and pressure on the brane. There is a contribution from higher order curvature 
corrections in these total quantities and hence in the modified Friedmann equation (23). In 
another words, the effect of the higher order curvature corrections in the brane action are 
interpreted as a contribution in the total energy density and pressure on the brane. On the 
other hand, the modified Raychaudhuri equation in this setup has the following form 



3m 



2 ' 



a 



+ p(--)-6r.(i/^ + A)] Uh-^) = -^(pM+p(-) + pM + p(--)) 



k 



X 



12mi 

k 



i2 



(24) 



2.2 Self-accelerating the normal branch with curvature correction 

In which follows we assume a spatially fiat FRW brane. In this case, we can rewrite the 
modified Friedmann equation (23) in the following form 

H^^^i^p(^)+p(GB)^^+e-. (25) 

where G = l/87rm| is the 4D gravitational constant. One of the simplest cosmological model 
that can exhibit acceleration of the universe is a de Sitter spacetime. So, it seems interesting 
to looking for such solutions in our model. With this motivation, we can probe the self- 
accelerating solution of the modified normal branch of this scenario. Assuming a de Sitter 
spacetime, i.e., looking for solutions of the form a{t) = aoCxp(iJo^), where Hq and Cq are 
constants, the modified Friedmann equation (25) in a de Sitter universe with Hubble rate 
Hq is given by 

+ Vo = constant. (26) 

where the subscript stands for quantities evaluated in the de Sitter phase and the scalar 
field is assumed to be constant in time. This relation shows that the matter energy density 
is constant by the virtue of the continuity equation (16), leads to a cosmological constant 
behavior for the matter sector (cj^™) = — 1) on the brane. In which follows, we neglect the 
matter content of the brane in our analysis of de Sitter brane. Using equation (23), the 
Hubble rate can be expressed as follows 



2rlH^ = 1 + + q/l + -^r^Vo. (27) 
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Note that the DGP braches can be recovered by imposing Vq = 0. Indeed, the case e = 1 
branch is corresponding to the de Sitter self- accelerating solution, and the e = —1 case 
corresponds to the normal branch leads to a Minkowski spacetime or the non-self-accelerating 
solution. In the presence of the curvature correction due to a Scalar-Gauss-Bonnet term in 
the action, the normal branch of the scenario self-accelerates too. It is important to note 
that the normal branch has not the difficulties of the self-accelerating one such as the ghost 
instabilities. In figure 1, we plot the behavior of the re-scaled squared Hubble parameter 
2rlHQ versus the re-scaled energy density ^^^^|^Vo for two branches of the solutions. This 
figure shows that the normal branch of the scenario exhibits a self-accelerating behavior 
because of the presence of an extra energy density originated from curvature corrections. 




3 f^n 

^ 3 

Figure 1: (a) Behavior of the re-scaled squared Hubble parameter 2r'^HQ versus the re-scaled en- 
ergy density ^^^^j^Vb for two branches of the scenario. The solid curve (with star) is corresponding 
to the normal branch and the dashed curve ( with circle) corresponds to the self-accelerating branch 
of the scenario, (b) A zoom on the normal branch of the DGP scenario with curvature effect. 



Up to this point, we have focused on the effects of a curvature correction term on the 
DGP model. Now we consider the effect of higher order embedding on the 4D Scalar-Gauss- 
Bonnet model. In other words, we study the effect of extra dimension on the cosmological 
dynamics on the brane. The Hubble parameter (27) can be expressed in a 4D regime as 
follows 

^0 = ah + • (28) 

where Hf^^ = ^Vq. This relat ion shows that the effect of an extra dimension is a shift on 
the 4D Hubble rate. In order to de Sitter brane to be close to the standard 4D regime ( with 
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i?(Q) ~ -f^(4)), the following condition should be satisfied 

< 1. (29) 



2.3 Phantom-Like Mimicry 

We can express the Priedmann equation (23) for the normal branch of this DGP-inspired 
scenario in a dimensionless form as follows 

= = + zf + noBil + ^)3(i+-''-) + 20,, 



-2^Qr.^Qm{l + zY + Qgb{1 + Zf(^+^GB) + (30) 

Q ^ , (m) Q ^ , (_GB) 

where = ^3^^ , Qgb = ^ sh'^ — ^f-c = 4;pF • Now we study the constraints 
imposed on this model. At redshift z = 0, equation (30) can be expressed as ( see Ref. [22] 
for a general framework in this direction) 

2 



^^nm + ^GB + ^r, - \l^r,j ■ (31) 

This expression indicates that Jl,, should be positive and Jl^ + f^cs + > too. Taking 
the square root of the last equation, there are two possibility for the sign of the quantity 
Q.jn + ^GB as follows 

+ ^GB > 0, (32) 

and 

-^r. <^m + ^GB < 0. (33) 

In the first case, taking the square root of the equation (31) yields 

1 + = ^Jn^ + nGB + ^r. ■ (34) 

Now by squaring both sides of this equation we find 

+ ^GB - 2^0^ = 1. (35) 

Similarly, for the second case the constraint equation is 

^^m + ^^GB + 2^= 1. (36) 

Note that the general relativistic limit can be recovered if we set fire — (o^ "^5 — 0)- 
this case equation (31) implies that flm + ^Igb = 1- 
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For a comparison between this DGP-inspired F{G, (f)) scenario and the well-known ADGP 
( see for instance [23]) and ACDM models, we study the expansion histories of these alterna- 
tive scenarios. We focus on the variation of the Luminosity Distance versus the redshift in 
these scenarios. The luminosity distance for a spatiaUy fiat universe is expressed as follows 



dL{z) r 

Jo 



dz 
Wz)' 



(37) 



where for the DGP-inspired F{G^ 0) scenario, the Hubble rate is given by equation (30). For 
ADGP model, H{z) is defined as [23] 



H{z) 



1 



roHo 



The luminosity distance of the ACDM model is given by 

Hgdz 



[1 + z 



I 

Jo 



flm{l + ZY + flf 



(38) 



(39) 



In figure (2) we show that the luminosity distance of the DGP-inspired scenario with cur- 
vature correction is closer to the ACDM model than the ADGP scenario. We note that 
this result is dependent on the choice of the model parameter space. Nevertheless, since the 
ACDM has a very good agreement with observational data, our analysis shows that with a 
suitable choice of the parameter space of the model, the DGP-inspired F{G, 0) scenario has 
better agreement with recent observation than ADGP. This result seems to be reasonable 
since the DGP-inspired F{G, 0) scenario has wider parameter space than the ADGP model. 



Now, we show that our model can lead to the phantom- like behavior on the brane without 

need to introduce any kind of the phantom scalar fields that violate the null energy condition. 
For this goal, the standard Friedman equation can be written as an effective form 



(40) 



where p^"*-* is the energy density of the standard matter and p^Jff^ is the energy density 
corresponding to dark energy. Comparison between the normal branch of equation (20) with 
this equation leads to the following relation for Pg//^ 



SnG (DE) _ (GB) , 3 
-^Pen -P +^ 



(41) 



where we have set SttG = 1. With the observable quantities, this equation is given by 

ijP^) . . 

= Qgb(1 + z)^(i+'^°«) + - 6^V^^(1 + zf + Qgb(1 + zf^^^-oB) + (42) 

tin 
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Figure 2: Luminosity distance versus the redshift for three alternative scenarios. In this figure we 
have set = 0.25, ilrc = 0.12, and ^Igb = —0.56. With these values, we calculate ^gb using the 
constraint equation (35). There is better agreement between DGP-inspired F{G,(j)) scenario and 
the ACDM at least in the parameter values adopted here. 

Note that pi^P and ujeff = yjj satisfy the continuity equation in the same way as the 
general relativity 

pi'^P + 3H{1 + ^e//)pi?f = 0. (43) 
In this respect, the effective equation of state of dark energy can be expressed as follows 

^^ff = -H , =- (44) 

pISP ^^m^ + + ^Gb(1 + Z^i^+^OB) + 

In figure 3 we plot the effective energy density and equation of state parameter of the model 
versus the redshift z. This figure shows a growing behavior of the effective energy density 
with decreasing of z; a requirement for phantom-like prescription. On the other hand, the 
effective equation of state parameter remains less than —1 in the phantom phase at future. 
Especially, the figure shows that the effective equation of state parameter has passed the 
phantom divide line uj^ff = — 1 in the near past. From figure 4 we can obtain appropriate 
ranges of z to fulfill the null energy condition, {pi^P ~^ -^efP — 0)- Clearly, this condition is 
satisfied in the present and future times at least in some subspaces of the model parameter 
space. As a result, the phantom-like prescription is realized in this DGP-inspired scenario 
without introduction of any phantom field neither in the bulk nor on the brane. 
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1.4- (a) (b) 




Figure 3: a) Variation of the effective dark energy density versus the redshift. The effective 
dark energy density grows by decreasing z and therefore shows a phantom-like behavior, b) The 
effective equation of state parameter versus the redshift. The braneworld universe has entered in 
the phantom phase in the near past and there is a smooth crossing of the phantom divide line in 
this setup. 




-0.8- 



Figure 4: The null energy condition is fulfilled for z < 1.1 . 
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3 An Extension with f{G) Term on the Brane 



In this section we consider a modified Gauss-Bonnet (modified-GB) term on the action of the 
DGP brane. The modified-GB gravity has been extensively studied in the hterature ( see for 
example [13] ). Especially in Ref. [24] it has been shown that this model may play the role 
of a gravitational alternative to dark energy proposal. On the other hand, the modificd-GB 
scenario has no ghosts and is stable, and it also passes the solar system tests. The action of 
the DGP-inspired modified Gauss-Bonnet gravity is given by 



Jbulk Z Jbrane 



-±R + mlK + f{G) + Cr^ 



(45) 



Note that the brane part of the action is dynamically equivalent to the Scalar-Gauss-Bonnet 
gravity with vanishing kinetic-energy term (see [24] for more details). In this action, f{G) 
is an arbitrary function of the Gauss-Bonnet invariant which is defined in equation (3). In 
the same fine as the previous section, the bulk-brane Einstein's equations calculated from 
action (45) are given by 

TIab - \gABn) + tuISa'^Sb'' {r,.. - \9^..R) S{y) = SA^'SB^T^Jiv), (46) 



where T^^, is the total energy- momentum on the brane defined as T^i, = T^™-* -|- T^p. The 
energy-momentum tensor Tj^^ is corresponding to the Gauss-Bonnet curvature on the brane 
and is given by 

rpi^i^if) ^ ]-gt'^f{G) - 2f{G)RR'"' + 4f{G)R^pR''f - 2f{G)R^""'^ R"" - 4f'{G)R^f""'Rp^ 

+2[V^VV'(G')]i?-25'^^[VV(G')]i?-4[VpV^/'(G')]i?'^^-4[VpVV'(G')]i?^''+4[VV(G')]i?''^ 

+4g''''[VpVJ'{G)]Rf"' - 4VpVJ'{G)]R''''''' , (47) 

where the prime denotes a derivative with respect to the G. From the equation (46), the 
field equations can be deduced as follows 

Gab = TIab - ^QabTI = 0, (48) 
and ^ ^ 

for the bulk and brane respectively. By imposing the junction conditions and using the FRW 
metric on the brane ( equations (17)- (19) ), we find the following Friedmann equation for 
the normal branch of the scenario 



VpM+p(/))-^. (50) 



3 ' y J- 
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Where we have set m| = 1. The energy density and pressure corresponding to the Gauss- 
Bonnet term can be defined as [22] 



and 



p(f) = Gf{G) - f{G) - 2A^h\2H^ + HH + AH^H)f"{G) (51) 

p(/) = f[G) - 24:'^H^(3H^ + 20H^H^ + 6H^ + AH^H + H^H)f{G) (52) 

respectively. Comparing equation (50) with the standard Priedmann equation, we can define 
an effective energy density as follows 

(53) 

where p^"*^ is the energy density of the standard matter and Pg//^ corresponds to the energy 
density of dark energy which is given by 

fT = P« - ^- (54) 

Now we assume f{G) to be defined as f{G) = fo\G\^ with constants /o and f3. For /9 < |, 
the f{G) term becomes dominant in the small curvature regime. We note that (3 can take 
essentially both positive and negative signs. Now the effective energy density can be rewritten 
in the following form 



& = hG^(hfi-l-{^y[2H' + HH + 'LH^H])-'-^. (55) 



As usual, the effective quantities piff ^-nd ujgff — yjj satisfy the continuity equation in 
the same way as the general relativity 

pi'^P + 3H{l + cu,ff)p^^P=0. (56) 

Considering a power-law expansion on the brane, we set the scale factor to be a{t) = a^t^^ 
with hf) = 1.2. This choice is reliable from physical grounds since it is corresponding to an 
accelerated expansion on the brane. Figure 5 shows the behavior of the effective dark energy 
density versus the redshift for two signs of /3. As this figure shows, although the effective 
energy density increases by decreasing z, its value remains positive only for negative values 
of f3. This property is illustrated explicitly in figure 5(6). In this figure we have plotted 
P^JfP versus (3 for various values of the redshift. The acceptable range of (3 can be deduced 

in a fascinating manner via imposing the null energy condition ( pi^P + PefP — 0) for this 
DGP- inspired modified-GB scenario as we have shown in figure 6. The null energy condition 
is fulfilled for a narrow strip defined as —0.25 < (3 < 0.05. 
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Figure 5: a) Variation of the effective energy density versus the redshift for (3 = —0.2 and (3 = 0.5. 
In both cases pfff^ increases by decreasing z, but its value always remains negative for positive 

values of /?. b) Variation of p^Jff^ versus f3 for z = ±0.25. For positive values of /? the effective 
dark energy is always negative and the phantom-like prescription breaks down in this case. 




■z=0.5 
. z= -0.5 



Figure 6: The null energy condition is fulfilled for a narrow strip —0.25 < /? < 0.05. 
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In figure 7(a), we have plotted the variation of the effective equation of state parameter 
versus the redshift for P = —0.2. Figure 7{b) shows that in the case of positive P, Uefj is in the 
quintessence phase (cUe// > — 1). For negative values of P, the effective equation of state lies 
in the phantom phase with {uJeff < — !)• A unified treatment of the above results extracted 
from figures shows that an effective phantom-like behavior can be realized in this setup in 
the region —0.25 < P < 0.05. Therefore, a DGP-inspired modified Gauss-Bonnet gravity 
has the potential to realize a phantom-like behavior in appropriate subspaces of the model 
parameter space. It is important to note that with scale factor defined as a = ao(^s ~ t)^°, 
an effective phantom phase can be obtained also with negative ho ( see for instance [24] 
and references therein). However, in our setup an effective phantom phase is realized with a 
positive value of ho by virtue of the induced gravity effects. 





Figure 7: a) Variation of the effective equation of state parameter versus the redshift for (3 = —0.2 
which remains in the phantom phase, b) Variation of the uJeff versus /? for z = ±0.25. In the case 
of positive /3, ujefj is larger than —1 but for negative values of /? the effective equation of state lies 
in the phantom phase with uJeff < — 1. 



4 Summary and Conclusion 

In this paper we have studied possible realization of the phantom-like behavior in a DGP- 
inspired scenario that curvature effects are taken into account by incorporation of the Gauss- 
Bonnet term in the brane part of the action. Although the Gauss-Bonnet invariant essentially 
should be considered in the bulk action, we can introduce it on the brane by considering a 
coupling between this invariant and a scalar field on the brane. In this viewpoint, one can 
consider possible modification of the induced gravity on the brane by incorporation of higher 
order curvature terms such as the Gauss-Bonnet type terms. Then, we have investigated the 
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phantom-like nature of this scenario with details. By the phantom-like behavior, we mean 
an effective energy density which is positive and grows with time and its equation of state 
parameter stays always less than —1, that is, the effective equation of state parameter is in 
the phantom region of the parameter space ( a;^^ < — 1 ). In this paper, firstly we have 
considered the normal branch of the DGP-inspired Scalar-Gauss-Bonnet gravity ( which is 
ghost free) and we have explored the cosmological implications of the model. Especially, we 
have shown that this branch, which is not self-accelerating in a pure DGP setup, can account 
for a self-accelerating behavior at the late times due to the Gauss-Bonnet curvature effect. 
As another important outcome of this setup, we have shown that an effective phantom-like 
behavior can be reahzed without need to introduce phantom fields neither on the brane nor 
in the bulk. The effective equation of state parameter of the model has a smooth crossing of 
the cosmological constant line at the near past and in the same way as observations indicate. 
In the last stage, we have considered a DGP-inspired model with modified-GB term in the 
brane action. In this case, the brane part of the action is dynamically equivalent to the 
Scalar-Gauss-Bonnet gravity with vanishing kinetic-energy term. We have shown that this 
scenario can also realize a phantom-like behavior by adopting cosmologically viable ansatz. 
The phantom-like prescription is realized in this setup for appropriate subspaces of the model 
parameter space ( for instance with —0.25 < /5 < 0.05). It is important to note that in the 
absence of the induced gravity effect, an effective phantom phase has been obtained with 
negative powers of the scale factor, that is, a phantom ansatz with ho < in a — ao(ts — 1)''°. 
Here we have shown that the presence of an induced gravity term on the brane leads to an 
effective phantom-like behavior with a positive values of ho which we call it a quintessence 
ansatz. Finally, we have shown that the luminosity distance of the DGP-inspired scenario 
with curvature correction is closer to the ACDM model than the ADGP scenario. However, 
this result is dependent on the choice of the model parameter space. Since the ACDM 
has very good agreement with observational data, our analysis shows that with a suitable 
choice of the parameter space of the model, the DGP-inspired F{G, 0) scenario has better 
agreement with recent observation than ADGP. This result seems to be reasonable since the 
DGP-inspired F{G, 4>) scenario has wider parameter space than the ADGP model. 
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